Measurable lattice effects on the charge and magnetic response in graphene 
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The simplest tight-binding model is used to study lattice effects on two properties of doped 
graphene: i) magnetic orbital susceptibility and ii) regular Friedel oscillations, both suppressed 
in the usual Dirac cone approximation, i) An exact expression for the tight-binding magnetic 
susceptibility is obtained, leading to orbital paramagnetism in graphene for a wide range of doping 
levels which is relevant when compared with other contributions, ii) Friedel oscillations in the coarse- 
grained charge response are considered numerically and analytically and an explicit expression for 
the response to lowest order in lattice effects is presented, showing the restoration of regular 2d 
behavior, but with strong sixfold anisotropy. 
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Introduction. The recent experimental realization [1] 
of graphene, the single layer honeycomb lattice of car- 
bon atoms that forms graphite, has unleashed an ex- 
plosion of activity. High expectations have been put on 
profiting from its peculiar electronic, mechanical, optical 
(and perhaps magnetic) properties, when tailored at the 
nanoscale [2] . The existence of linearly dispersing bands 
around two nodal points [massless Dirac fermions with 
velocity v ^ 10^ m/s) form the basis of graphene's most 
notable electronic properties [3]. 

Many theoretical studies of graphene are done within 
scaling limit or Dirac cone approximation, that is, as- 
suming strictly linear energy dispersion around the nodal 
points. Although this approach is successful in explain- 
ing many experimental facts, it has limitations too. Ob- 
vious examples are provided by magnitudes for which the 
Dirac cone approximation provides a null result. In this 
paper we are concerned with two such magnitudes: i) 
the magnetic susceptibility and ii) regular Friedel oscil- 
lations, both rendered zero at finite doping in the scaling 
limit. 

i) Strong and peculiar diamagnctism characterizes 
graphene, as first discussed by McClure to explain 
graphite. Nature's best diamagnet. He found that, for 
the two-dimensional Dirac model, the diamagnetic sus- 
ceptibility was given by a delta-function of the chemical 
potential [4]. This result implies that there is no mag- 
netic response when the chemical potential is shifted from 
the neutrality point. This is in clear contrast with recent 
experimental findings of paramagnetism in graphene [5] . 

Here we will show that lattice effects, neglected in 
the scaling limit, render finite and sizable the magnetic 
response. We employ the formalism of Fukuyama [6], 
whose original formula was first applied to graphite [7] 
and subsequently to graphite intercalated compounds 
[8, 9], which is here extended by an additional term 
required to provide the exact susceptibility for a gen- 
eral tight-binding model. The magnetic response for 
arbitrary chemical potential is obtained, finding orbital 



paramagnetism (OP) over a wide range of fillings. Its 
value close to the neutrality point is compared with other 
sources (core diamagnctism, spin paramagnetism and in- 
teraction's induced orbital paramagnetism) and shown to 
be a relevant contribution. 

ii) Graphene's charge response around localized pertur- 
bations is also peculiar [10-12]. While ordinary 2d sys- 
tems show the familiar Friedel 2kp oscillations decaying 
as I/t"^, graphene's coarse-^ramerf response in the scaling 
limit does so but with an additional power. Graphene's 
lack of regular 2d Friedel oscillations is linked to isospin 
(or chiral) conservation and thus provides the possibility 
of direct observation of the nature of graphene's excita- 
tions [13]. 

Here we show, numerically and analytically, that lat- 
tice effects restore the standard 2d behavior. An explicit 
expression for the charge response is obtained to lowest 
order in lattice effects, exhibiting the usual l/r^ decay 
and a pronounced sixfold anisotropy, with maxima along 
the bond's directions. 

One might ask why we treat two at first sight such 
distinct topics on the same footing. The reason is that 
within the Dirac cone approximation the static transverse 
current-current as well as the charge-charge correlation 
function yield xio) = a + bq^ + ... with b ^ 5{Ep) and Ep 
the Fermi level. Lattice contributions to the response, 
given by ^ for finite filling factor, are thus 

suppressed in the same peculiar way. 

Tight-binding model and Dirac cone approximation. 
We describe graphene by the simplest tight-binding 



H.c. with hop- 



Hamiltonian H 
ping amplitude —t between nearest neighbor atoms in 
A and B sublattices joined by vectors = (0, a) and its 
±120°-rotated versions 62,3- The spectrum is E'^{k) — 
±|t5'(fe)|, where S{k) = ^ ■ exp(zfe • 8i). It develops 
a well-known linear dispersion E^{k) — ±hvk' with 
k' = k Ki 2 around two points in the Brillouin zone, 

For the orbital magnetic response, the Dirac cone ap- 
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proximation leads to a diamagnetic susceptibility de- 
pending on Fermi level Ep as [4] 
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(1) 



with vacuum permeability fig (SI units), spin and valley 
degeneracies Ss = 9v = 2, and unit charge e. The pecu- 
liar relation, known to be at the basis of graphite's large 
diamagnetism [7] , implies that doped graphene shows no 
magnetic orbital response at finite doping. 

Also the charge response to a local impurity V = uS{r) 
shows peculiar behavior since the first (/-derivative of its 
susceptibility is continuous at q = 2k p. This results in 
an anomalous decay of the Friedel oscillations which in 
terms of the electronic carrier density reads [10] 

Sp{r) upe cos(2fci?r) 



kpr > 1 . 



(2) 
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In what follows, we will show that both results are sub- 
stantially altered when lattice effects are included. 

Magnetic response. Fukuyama [6] has provided a con- 
venient expression for the orbital magnetic susceptibility 
in non-interacting systems. The formula is exact only for 
a Hamiltonian of the canonical form H = ^ — \-V(R). 
Here we adapt Fukuyama's procedure to obtain the exact 
orbital response of a tight-binding Hamiltonian. Employ- 
ing the current operator of the tight-binding model given 
in Ref. [14], linear response theory yields the following 
expression for the orbital susceptibility [15]: 
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with the 2x2 matrix Ok (E) = (E 
'^kHk, and Hk given by 

-t S{k)\ 
-t iSik)r ) 

This expression for Xorb transcends the model of the 
initial Hamiltonian, and turns out to be correct for any 
tight-binding system. Eq. (3) coincides with Fukuyama's 
original result [6] except for the second term. The differ- 
ence stems from the standard isotropic P-dependence in 
Fukuyama's H — ^ + ^(-R)i where one has |^ = 0. 
Such cancellation does not apply in a generic tight- 
binding case. 

The above formula has been applied to numerically 
calculate the orbital susceptibility in graphene as a func- 
tion of Fermi level Ep- The most prominent feature is, of 
course, the delta function at the band center of Eq. (1), 
which comes with the known analytical value. We can 
thus extract the lattice contribution from the numerical 
results by writing 

9s9v 2 2; 
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FIG. 1. Continuous line: lattice contribution to the orbital 
magnetic susceptibility (Axort) in units of Xo = Hofi^^e'^\t\a'^ ■ 
Dashed hne: Pauh's spin paramagnetic contribution. Dashed- 
dotted hne: orbital magnetic susceptibility from Coulomb 
electron-electron interactions (Eq. (6)). The arrow marks 
an estimate of the absolute value of the diamagnetic contri- 
bution from core electrons [17]. 



The calculated lattice contribution A.Xorb is plotted in 
Fig. 1 in units of Xo = /io^~^e^|t|a^. The lattice origin 
of Axorb becomes evident if one artificially sets a — 
and t oo while at cc v — constant (scaling limit): then 
Axorb 0, leaving the Dirac cone result (Eq. (1)) as 
the sole response. 

Discussion. The orbital response is usually associated 
with diamagnetism, so a noteworthy aspect of the lattice 
contribution in Fig. 1 is its paramagnetic character over 
most of the band, even diverging at the van Hove points 
[16]. From Eq. (3), one finds the following sum rule: 
JdEp Xorb{Ep) = 0. The existence of OP is thus a 
necessary consequence to cancel the large diamagnetic 
contribution of the scaling limit at the band center, see 
Eq. (5). Only at the band edges Landau diamagnetism 
emerges, as expected. 

Now, we compare lattice's OP with other contribu- 
tions in the region Ep ~ 0, relevant for gate-doped 
graphene. Within our non-interacting model, the only 
remaining magnetic contribution is Pauli's spin param- 
agnetism, given by Xspm/tJ^o = '^fJ'sPi^p) where pB is 
the Bohr magneton and p{Ep) is the density of states 
(per spin). This spin contribution is plotted in Fig. 1, 
where it is seen that it cannot compete with the domi- 
nant orbital contribution for low carrier densities pe- 

Core electrons, not considered in our ppvr band Hamil- 
tonian, are another source of (dia)magnetic response. 
The estimate of Ref. [17], Xcore ^ —4.8 x 10~^emu/mol, 
translates into Xcore ^ — 0.13xo, a value marked with 
an arrow in the scale of Fig. 1. Again, the orbital con- 
tribution for low Pe is comparable or greater than this 
estimate of core diamagnetism. 

Electron-electron interaction is the major ingredient 
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left out in our model. Recently, its effect on the mag- 
netic response has been calculated to first order within 
the Dirac cone approximation [18]. Two scenarios have 
been considered: Thomas-Fermi screened Coulomb in- 
teraction and contact (Hubbard-like) interaction, both 
leading to OP at finite doping. For the Coulomb case, 
the interaction's contribution to the susceptibility can be 
written as [18] 



Xe 
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Xo 



Ep/\t\ 



(6) 



with an interaction dependent constant C ~ 0.04, suit- 
able for graphene over Si02. Xee is compared with the 
lattice contribution in Fig. 1. Graphene's poor screen- 
ing causes the divergence of Xee when Ep — )■ 0, making 
this contribution dominant. But even in this unfavorable 
case, the lattice contribution is not negligible. For in- 
stance, Axorb ^ Xee for doping levels n> 3.7 X 10^'^ cm~^, 
and for a typical doping n 4 x 10^^ cm~^ one has 

Xee ^ SAXorb- 

Finally, if (by some external means) screening were 
truly effective so that interactions could be described by 
a contact term Vo, Rcf. [18] provides the following ex- 
pression for the interaction's promoted orbital paramag- 
netism: 
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where we have written the interaction as — 
with a Hubbard-like energy U, and area per unit cell 
Ac = ■ In order to compare this contribution with 

that of the lattice, wo calculate the size of U required for 
Xee of Eq. (7) to match the lattice contribution around 
Ep ~ 0. The answer turns out to be 7\t\, a value 
substantially larger than current estimates for graplicnc. 
This implies that, in any reasonable scenario of contact 
interactions, the lattice orbital contribution to paramag- 
netism would be the key player in the magnetic response 
of doped graphene. 

Let us close with a remark on the effect of next-nearest 
neighbor hopping, temperature and disorder. As was al- 
ready noted in Ref. [9], t' « O.lt leads to a consider- 
able electron-hole asymmetry in the magnetic response; 
the above qualitative discussion on the relevance of the 
several contributions, though, is not altered. Temper- 
ature and disorder [19] broadens the diamagnetic delta- 
peak, such that lattice effects gradually lose relevance at 
a given temperature or disorder when the chemical po- 
tential decreases to zero. 

Charge Response. The linear, static charge response of 
graphene is given by 

mr.^ 1 f (I. I. , np{E^{k,))-np{E^'{k2)) 

(8) 



with k\=k2-\- and the prefactor f± 

/±(fer,fc2) = -±-Re^^^^^ 



(9) 



Friedel oscillations arc caused by intraband transitions 
(-h- sign in Eq. (9)) with q/2 ~ kp, the Fermi wavevec- 
tor (measured from the Dirac point). To understand 
graphene's peculiarity, let us set /+ = 1 in Eq. (8), 
and call the associated Lindhard-like response 11. Then, 
the dominant singularity in H, corresponding to transi- 
tions across the Fermi surface, leads to the prototypical 
2d square-root behavior 
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where G is Heaviside's function and we have ignored any 
distortion of the isotropic Fermi surface around the two 
Dirac points. Within the Dirac cone approximation, the 
prefactor /+ in Eq. (10) crucially vanishes for states on 
opposite sides of the Fermi surface. The physical inter- 
pretation is well known [10]: the involved states have 
opposite chirality and cannot be coupled by a perturba- 
tion diagonal in sublattice index. Nevertheless, this ex- 
act cancellation of /+ holds true only in the scaling limit 
kpa — 0, and a finite value of kpa renders /+ finite, 
something we generically label as lattice effect. 

To see if this square root behavior is present also for 
the true prefactor /+ as given in Eq. (8), we numerically 
analyze the response derivatives which we conveniently 
write as 



dUiq) 



-J^Im jdEnp{E)Y, (11) 



T,{Qk{E)^ZQk{E){gk-q{E) - gk+q{E))} . 

We observe a clear anisotropy with a pronoimced spike 
that hints at a singular behavior for the results in the y 
direction, absent in the x direction, where the behavior is 
closer to that of the (analytically known) Dirac cone ap- 
proximation [15]. The numerical results strongly suggest 
the restoration of a regular 2d response but with strong 
anisotropy. This is confirmed by the analytical treatment 
that follows. 

Now we obtain analytically the charge response to low- 
est order in lattice effects. We start with the determi- 
nation of the prefactor /+(q), that is, Eq. (9) for two 
states on opposite sides of the Fermi surface: k2 and 
fci = k2 + q, such that the vector q corresponds to the 
square-root singularity in the response. The latter condi- 
tion requiring the g-linked portions of the Fermi surface 
to be parallel. Upon a Jacobi-Anger expansion of the 
terms exp(ife • Si), the structure factor can be written as 
S{k) = 3E„ J_i+3„(fc'a)e'(-i+3")'^ where J„ are Bessel 
functions of the first kind, and the separation from the 
Dirac point is k' = k — Ki, with polar coordinates (fc', (p). 
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To lowest order in lattice effects, only J_i = — Ji and J2 
are to be retained. Then, the requirement = const 

leads to the following expression for the Fermi surface in 
polar coordinates: 



k'p{(j)) = kf 



1 , ^^"^ 

1 H — — cos(30j 



Oikpaf 



(12) 



where we have parametrized the Fermi energy by the 
would-be Fermi wave vector in the isotropic limit: Ep = 
tivkp- 

To lowest order, the condition of parallel pieces of 
Fermi surface leads to the following relation between po- 
lar angles {(f)i,2) of the involved fe-points: (pi = (pq — 6<p 
and (/)2 = 0q + TT + (50 with the lattice correction Sep = 
I {kpa) sm{3(j)q) and (pq the polar angle of the joining 
vector q with modulus q = 2A;ir(l + 0{kpaf'). We can 
now write the phase of S{k') as 



S{k') 
\S{k')\ 



_^-i{'t>+e) 



(13) 



where cp is the polar angle of k' and 9 is the lattice 
correction to that phase given to lowest order hy 9 = 
j{kFa) sin(30). This leads to the final expression for 
the prefactor 



{kpaf 



(1 - cos{%(pq)) + 0{kFaf 



(14) 



where {2kp, (pq) arc the polar coordinates of q and which 
holds for both Dirac points. We note that the result 
of Eq. (14), although the lowest finite order m a, kp 
expansion, already represents an excellent approximation 
for sizable Fermi levels well within the range of gate- 
voltage doped graphene's samples [15]. 

Combining Eqs. (10) and (14), the dominant singular- 
ity of the true response is 
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We can now determine the density response associated to 
a local perturbing potential V = u5{r) given by = 
■^2^ / d^q e"^ 'T[(q'). The remaining integral is obtained 
from standard techniques [20], leading to the following 
asymptotic behavior for Priedel oscillations: 



5p{r) {kpaf 



Pe 
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(1 — C0s{6(pr)) 



upe sm{2kpr) 
Ep {2kpry 



(16) 



While its r^^ behavior is standard for a 2d system, its 
true origin as a lattice contribution to an otherwise null 
result (to order r~^) is revealed by its amplitude, van- 
ishing as {kpa)^ in the scaling limit a — >■ 0, and by its 
anisotropy reflecting the sixfold symmetry of the lattice. 

Comparing Eq. (16) with the result coming from the 
Dirac cone approximation, Eq. (2), we first notice the 



phase shift of tt/2. We further find for the crossover 
length between anomalous and regular Friedel oscilla- 
tions ~ kp^{kpa)~^. For kp = lnm~^, we have 
Tc ^ lOOnm which corresponds to an impurity concen- 
tration of rii ~ 10^°cm~-^, recently found to be the in- 
trinsic concentration of local impurities in graphene [21]. 
We thus expect Friedel oscillations to show anisotropic 
behavior and modify the RKKY-interactions for doping 
levels Ep>0.5eV. 

Summary. The simplest tight-binding model has been 
employed to study the lattice contribution to the mag- 
netic susceptibility and (coarse-grained) charge response 
of doped graphene, for which the Dirac cone approxima- 
tion produces a null result. The lattice magnetic response 
shows orbital paramagnetism for a wide range of filling 
factors, representing a relevant contribution when com- 
pared to other sources such as core diamagnetism, spin 
paramagnetism and electron-electron interaction induced 
orbital paramagnetism. Lattice effects restore the 2d 
regular behavior for the coarse-grained charge response, 
with Friedel oscillations decaying as r~^ but with pro- 
nounced sixfold anisotropy, with maxima along bond's 
directions. For clean samples with impurity concentra- 
tions of rii 10^°cm~^ they become relevant for doping 
levels Ep>0.5eV. 
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